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ABSTRACT 

The matrix element of the kinetic energy operator between B meson states 
is computed by means of a QCD relativistic potential model, with the result: 

= 0.66 GeV 2 . A comparison with the outcome of other theoretical ap- 
proaches and a discussion of the phenomenological implications of this result 
are carried out. 
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1 Introduction 



In the last two decades, the study of hadronic processes involving heavy quarks has 
attracted continuous interest both from experimental and theoretical sides. The main 
theoretical achievements have been obtained in the framework of Heavy Quark Effective 
Theory (HQET) JTJ, which describes the dynamics of heavy hadrons, i.e. hadrons contain- 
ing a heavy quark Q, when mq — > oo. The theory is based upon an effective lagrangian 
written in terms of effective fields, which is a systematic expansion in the inverse powers 
of the heavy quark mass mq. In particular, it has been pointed out that the expansion 
in the inverse powers of mq is nothing else but an application of Operator Product Ex- 
pansion (OPE) in the sector of heavy hadrons ||. The leading order effective Lagrangian 
displays heavy quark spin and flavour symmetries which are not present in full QCD. 
These symmetries are no longer conserved at the next-to-leading order, and the 0(^^j 
lagrangian reads as follows: 

C = h v iv- Dh v + - h v [D 2 - (v ■ Df]h v + -^±-h v ^—h v + O(-) . (1) 

Zrriq ZfTlq Z m Q 

In the mq — > oo limit the field h v is related to the heavy quark field Q by: Q = 
e~ imQV ' x h v (x), where is the heavy quark four-velocity which, in the mq — > oo limit, 
coincides with the one of the hadron 0. In the hadron rest frame, the first of the two 
next-to-leading order operators appearing in (HD is the heavy quark non relativistic ki- 
netic energy due to its residual motion, while the second one is the Pauli chromomagnetic 
interaction operator; they correspond in the Wilson expansion to dimension 5 operators. 
Their matrix elements can be parametrized as follows: 

„2,„ n < H Q \h v {iV) 2 h v \H Q > 

^(H Q ) = _ (2) 

»g{H q ) = — (3) 

where Hq denotes generically a hadron containing the heavy quark Q and = d^ — igA^ 
is the covariant derivative. The normalization: < Hq\Hq >= 2Mh q is understood. 

Since Hq represents the chromomagnetic interaction between the heavy quark spin 
sq and the light cloud total angular momentum sg, it can be obtained from the mea- 
sured hyperfme mass splitting, when available. Its general expression reads: ij,q(Hq) = 



1 



— 2[J(J + 1) — |]A2, where J is the total spin and A2 is independent of the heavy quark 
mass. Therefore, in case of B mesons, we have: 

nUM h ) = d -f{Ml*-Ml). (4) 

where M b = B,B* and du = 3 in the pseudoscalar case, du = — 1 in the vector case 
(hence, from experimental data [|: n%(B) ~ 0.36 GeV 2 , n%(B*) ~ -0.12 GeV 2 ). On the 
other hand, it is expected to be zero for all baryons whose light cloud is in a se = state, 
such as Aq, Sq, while it should not vanish in the case of Qq for which se = 1, though it 
is not experimentally known, yet. Moreover, the mass splitting has been measured in the 
case of E b : M s * - M Sb = 56 ± 16 MeV §. 

H 2 , represents the average square momentum carried by the heavy quark inside the 
hadron, that is, modulo a factor 2m<g, its non relativistic kinetic energy. 

These quantities are interesting for several reasons. Heavy hadrons masses are ex- 
pected to scale with rng as: 

M Ho=mc + A + 4^ + ... (5) 

A represents the difference between the mass of the hadron and that of the heavy quark 
in the mg — > 00 limit. In this limit, it can be related to the trace anomaly of QCD ||: 
A = < H Q \^-G^G^\H Q >, where p is the Gell-Mann-Low function. Moreover, 

if the inclusive semileptonic width of a heavy hadron is calculated by an expansion in the 
powers of the following results are found: the leading term of the expansion coincides 
with the free quark decay rate (spectator model); no corrections of order ^ affect the 
rate; the ^ corrections depend on fi 2 and /j/q 0. The absence of — corrections is 
sometimes referred to as CGG/BUV theorem. As a consequence, these paramenters enter 
in the ratio of hadron lifetimes and in the lepton spectrum in inclusive transitions, which 
in principle are quantities directely comparable with experimental data []. 

It is worth noticing that fi 2 and \j? g are the matrix elements of operators which are 
not sensitive to the light quarks flavour, i.e. they are SU(3) singlet operators. The SU (3) 
breaking effects emerge at —3- level, due to four-quark operators; their matrix elements 

m Q 

can be estimated by factorization in the case of mesons, and, in the case of baryons, by 
constituent quark models || or field theoretical approaches, for example QCD sum rules 



2 A critical analysis of such a procedure can be found in 
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In this work we calculate [i 2 in the case of B mesons by means of a relativistic potential 
model. In the next section, after describing the relevant features of the model, we will 
present our results. Phenomenological implications and comparison with other approaches 
will be provided in section 3. Finally, we will draw our conclusions. 



2 Method and Results 



We will describe in the following the relativistic quark model used to compute [i 2 . Within 
this model, the state of a pseudoscalar (bq a ) meson is written in terms of a wave function 
ipB and of quark and antiquark creation operators; in the meson rest frame it reads: 



B a >= J dk Mk) b\k,r,a) 4(-fc,a,/3)|0 > 



(6) 



In (Q) a, (3 are colour indices, r,s are spin indices and a is a light flavour index; the 
operator creates the b quark with momentum: k, while d\ creates a q a antiquark with 
momentum —k. The wave function ipg is obtained as a solution of a Salpeter equation 
n | which takes into account relativistic effects in the quark kinematics: 



{\jk 2 +ml + ^/k 2 + m 2 a - M B }4j B {k) + J dk'V{k, k')tfj B {k') = 



(7) 



Eq. (0) stems from the quark-antiquark Bethe-Salpeter equation in the approximation 
of istantaneous interaction. The interquark potential V is represented by the Richardson 
potential |T3], which reads in the r— space: 

8vr . r . /(Ar) 



V(r) 



-A 



33 - 2n/ 

A is a parameter (chosen at the value A = 
flavours, and the function f(t) is given by: 

sin(qt) 



Ar 



Ar 

397 MeV M 



71/ is the number of active 



/(*) 



4 
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dq- 



q 



1 



ln(l + q 2 ) q 2 



(9) 



The potential (^) is linear at large distances in order to assure QCD confinement; at short 
distances it behaves as — with a s (r) logarithmically decreasing with the distance 
r to reproduce the asymptotic freedom property of QCD. Spin interaction effects are 
neglected since in the case of heavy mesons the chromomagnetic coupling is of order tbq 1 . 
The masses of the constituent quarks are fixed in such a way that the meson spectrum 



of the charmonium, the bottomomium and of the heavy-light systems is reproduced: the 
fitted values are: m& = 4.89 GeV, m q = m u = ma = 0.038 GeV. The mass of the charm 
quark has also been fixed: m c = 1.452 GeV. Finally, the wave function ip B is covariantly 
normalized: 

j^— J dk\Mk)\ 2 = 2M B . (10) 
In the B meson rest frame it is useful to define also the reduced wave function u B (k) 

(k = \k\y. 

k^ B (k) 

u B {k) = -/=- (11) 

which is normalized as: dk\u B (k)\ 2 = 2M B . The function u B (k) can be obtained by 
numerically solving eq. (^) using the Multhopp method described in [14, [13|]. The B— 



meson wave function, together with the wave functions computed for the other mesonic 
states analyzed within this framework, is the main outcome of the model; by using it a 
number of hadronic quantities characterizing the B system have been computed, such as 
semileptonic form factors, leptonic decay constants and strong coupling constants [JT3J, [15 



By writing the heavy field, in the expression of the kinetic energy operator, in terms 
of creation operators, and by exploiting usual anticommutation relation and the normal- 
ization condition for the wave function, we obtain a simple expression for the 6-quark 
average momentum squared in the B meson: 

K imt) = y . , 12) 

J dk \u B {k)\ 2 

The parameter fi 2 (Bd) should coincide with K(rrib) in the limit m& — > oo. In order to 
perform such limit, we compute K{m^) for several values of m^, using the wave function 
obtained in corrispondence to the appropriate value of m&; the result is plotted in Fig. 1. 
We then extrapolate the resulting curve to mj — > oo according to the expression: 

tf(m 6 ) = A 4(£ d )(l + — ). (13) 

The result is: 

fil(B d ) = 0.66 GeV 2 (14) 

a = -1.54 GeV. (15) 

It should be noticed that this extrapolation procedure cannot be avoided in our approach 
since the resolution of the Salpeter equation can only be performed by numerical methods. 
Let us observe that the result: a = —1.54 GeV suggests that for the charm the — 



corrections play still an important role. On the other hand, the result obtained using the 
value of rrib fixed within our model, i.e. mj, = 4.89 GeV, is: K = 0.46 GeV 2 , which means 
that for the b quark the finite mass result differs from the asymptotic one at the level of 
30%. Moreover, eq. ([L4]) agrees quite well with the QCD sum rule result |2£| . 



This result can be translated in a determination of the parameter A. As a matter of fact, 
exploiting eqs. (f|) and (H) the following relation can be derived: 

3M B * +M B T ul 

M B = ?- *=m b + A + -t^ 16 

4 2m h 



which, using the value in eq.(|T4]), gives: 

A = 0.35 GeV . (17) 

The results (14) and ( |17|) are also consistent with (||) if one neglects the higher order 
corrections in l/m& and uses the experimental value of fiQ^Bj). 



Let us briefly discuss the uncertainties of the result in (|14T) . They mainly depend on 
the computed wave function Us{k), and, therefore, they can be estimated by modifying 
the shape of the wave function. We perform this analysis by comparing the outcome in 



(14) with the result of a similar constituent quark model, such as the one in ref. [16] 



where the value of /j, 2 (Bd) is obtained by two independent methods. The first one consists 
in using the Altarelli et al. (ACCMM) model [17], where the heavy quark momentum 
distribution inside the B meson is assumed to be gaussian: 4>{\k\, P F ) = exp(^ — , 
P F being a parameter. In terms of Pp the kinetic energy fil(B d ) reads: /i^(£>d) = l-Pf. 
The value of Pp has been obtained in [16| by a comparison with recent CLEO data [Tj| 



on the inclusive B — > Xiv semileptonic rate: Pp = 0.54 ±0^5 GeV, which corresponds 
to: fi 2 ~ 0.44 GeV 2 . The second method consists in using a quark model, which gives 
Pp = 0.5 — 0.6 GeV, i. e. /i^ = 0.375 — 0.54 GeV 2 . We may compare our result with 
the range of values quoted above: /i^ = 0.375 — 0.54 GeV, observing that the comparison 
must be performed with our finite mass result, since the wave function used in [jl6l has 



been obtained for real values of rrib. By this comparison, we can conservatively conclude 
that our result is affected by an error of 20% related to the shape of the B meson wave 
function. 
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3 Phenomenological implications 



Various determinations of the value of fi 2 (Bd) exist in the literature they are collected 
in Table I. 



The analyses in refs. ||16|| , P(|-||25|| consist in various attempts to extract or to put 



constraints on nl(B d ) from experimental data. In refs. PU| , f23| experimental data on 
semileptonic B and D decays are compared to theoretical predictions to extract fjL^ as a 
function of A ; in particular, in ref. |^U| the QCD sum rule result [EBJ: A = 570 ± 70 MeV 
is used to constrain /i 2 in the range: 0.1 < fi 2 < 1.5 GeV 2 . A similar approach is employed 
in ref. where it is stressed the possibility of obtaining A and fi 2 from the moments 
of the photon spectrum in the decay B —>■ X s j. 

QCD sum rules have been applied in refs. |27j and |2£| to determine fi 2 . Moreover, a great 



deal of works have been devoted to further constrain theoretically [i 2 . A field theoretical 
approach has been applied in ref. || (confirmed by a quantum mechanical approach in 
p9|, p0|1 ) to state the inequality: /i^ > p? G []. Moreover, in a theoretical argument has 



been given to confirm and strengthen the bound, giving fi 2 > 0.45 GeV 2 in the case of 
B mesons. This argument is based upon the possibility of extracting /i 2 from the slope 
of the Isgur-Wise function which is related, at the leading 1/rriQ order, to the differential 
decay rate: ^{B — > D*lv) which can be obtained from experimental data . 
Finally, very recently /i 2 has been computed on the lattice |33| . 



On one hand, this variety of results suggests that further theoretical analyses of /j, 2 are 
required, on the other it shows that the experimental determinations are a hard task. The 
main difficulty lies in the smallness of the parameter /ijj. and in the fact that it appears in 
physically measurable quantities always in connection with quantities that are determined 
in more or less broad range of values, such as A and the quark masses m c , mi,. 

As an example, we may consider the role played by fi 2 in the B semileptonic branching 
ratio, a well known problem in B physics, since theoretical estimates are still larger than 
experimental data. The most recent experimental measurement has been performed by 



CLEO Collaboration || giving: B(B -> Xlv) = 10.49 ± 0.17 ± 0.43 %. 



From the point of view of the — expansion, the general procedure to determine an 
inclusive quantity consists in applying the OPE to the forward matrix element of a weak 



3 In the baryon sector, a result for the has been obtained by QCD sum rules: /z^.(A&) ~ 0.6 GeV 2 
[ill with an estimated uncertainty of 30%. 

4 In |2^| the bound has been criticized, stating that the role of radiative corrections prevents any 
possibility of constraing /i^.. 
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transition operator . The resulting expression for the lepton spectrum in B semileptonic 
decay, derived in reads: 

dV 

dy 



T 9(1 -y- p)2y 2 { (1 - ff{l + 2/)(2 - y) + (1 - /) 3 (1 - y) 



+(!• 
-5 



/)[(!-/) (2 + 



1 - /) 2 (1 + 2f)y + f -{l - /)(10y - 8y 2 ) + ^(3 - 4/)(2y 2 - y 3 )] ^} 



2/ + ^/y) - -[2y + /(12 - 12y + 5y 2 )] 
3 1 P 

f 4 



Pg_ 

Pi 



where: 



G\m\ 

1927T 3 



w. 



qb\ 



f 



ni 



2 ■ 



(19) 



and m q is the mass of the final quark q f\. 



Let us consider the case B — > X c tv. Eq. (^) depends on the value of the charm quark 
mass m c , which is fixed in the potential model to the value m c = 1.452 GeV by fitting the 
whole charmonium spectrum . Using V& = 0.04 and our result (14) we obtain the lepton 
spectrum displayed in fig. 2. This curve must be compared to the experimental distibution 
Since on the theoretical side, heavy quark mass expansion is unreliable for large 



in 
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lepton energies Et>2 GeV [35], and, in the experimental analyses, high energy leptons 
must be selected in order to subtract the background of secondaries, the comparison 
theory-experiment can be performed only in a selected window of lepton energies. For 
example, as it has been done in refs. 



Ri 



E, >1.5 GeV E i-dW e dEj 



25] , one can use the ratios: 
I 



dT 



dE, 



>E t >l.5GeV dE L 

where the dependence on the overall factor |Vd,| 2 m^ cancels 



E t >1.7GeV dE e dE £ 
JEi>1.5GeV dE~ f dEl 



(20) 



Using the experimental 

results: R\ = 1.7830 R 2 = 0.6108 for the ratios in ( |2"0"D the values of p 2 displayed in 

Table I have been obtained [24], [25] 0. Using the formulae for Ri^ in [Q and our results 

flilif) , (|l7|), we obtain: R\ = 1.733 R 2 = 0.559. However, as already pointed out in fl25i , 

the parameters A, /x 2 enter in i? 12 as power corrections and represent a small effect in 

(20), so that very small changes in the theoretical or in the experimental expressions for 

Rip would shift the values of A, /i 2 towards very different results. One of such changes 

could be related to a different estimate of the secondary electron background, or, from the 

5 The expression of 0(a s ) corrections to cq. ( [l8|) can be found in |3q| , 
6 The formulae in |24 for R12 include also 0(a s ) corrections. 



theoretical side, to the next order perturbative corrections, whose size is difficult to assess. 
Therefore, one cannot avoid to conclude that the accuracy of a single determination from 
experimental data is difficult to check. A set of independent measurements, from different 
channels, should be used, and an accurate cross check of the errors should be performed 
to detect the value of the parameter /i^. 

Acknowledgments 

I thank P. Colangelo and G. Nardulli for interesting discussions. 



8 



References 

[1] For a recent review see for example: B.Grinstein, An Introduction to Heavy Mesons, 
published in: Mexican School 1994, pag. 122 and references therein. 

[2] For a comprehensive review on this subject see for example: M.A. Shifman, Lectures 
on Heavy Quarks in Quantum chromodynamics, preprint TPI-MINN-95/31-T, UMN- 
TH-1413-95, |hep-ph/9510377| . 

[3] H. Georgi, Phys. Lett. B240 447 (1990). 

[4] Review of particle properties, Phys. Rev. D50 1173 (1994). 

[5] P. Abreu et al, DELPHI Collab., preprint DELPHI 95-107 PHYS 542 (1995). 

[6] I. Bigi, M. A. Shifman, N. Uraltsev and A. Vainshtein, Phys. Rev. D52 196 (1995). 

[7] J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B247 399 (1990); I. Bigi, N. Uraltsev 
and A. Vainshtein, Phys. Lett. B293 430 (1992); (E) B297 477 (1993). 

G. Altarelli, G. Martinelli, S. Petrarca and F. Rapuano, preprint CERN-TH/96-77 , 
|hep-ph/ 9604202 . 



[9] N.Bilic, B.Guberina and J.Trampetic, Nucl. Phys. B248 261 (1984). 

[10] P. Colangelo and F. De Fazio, preprint BARI-TH/96-230, |hep-ph/ 9604425 

[11] E. E. Salpeter, Phys. Rev. 87 328 (1952). 

[12] J. L. Richardson, Phys. Lett. B82 272 (1979). 

[13] P. Colangelo, G. Nardulli and M. Pietroni, Phys. Rev. D43 3002 (1991). 

[14] K. Karamcheti, Principles of ideal fluid aerodynamics (Wiley, New York, 1966). 

[15] P. Cea, P. Colangelo, L. Cosmai and G. Nardulli, Phys. Lett. B206 691 (1988). 

[16] D. S. Hwang, C. S. Kim and W. Namgung, preprint KEK-TH-473, |hep-ph/ 96042251 . 

[17] G. Altarelli, N. Cabibbo, G. Corbo, L. Maiani and G. Martinelli, Nucl. Phys. B208 
365 (1982). 

[18] B. Barish et al, CLEO Collab., Phys. Rev. Lett. 76 1570 (1996). 

9 



[19] P. Colangelo, C. A. Dominguez, G. Nardulli and N. Paver, Phys. Rev. D54 4622 
(1996). 



[20 
[21 
[22 
[23 
[24 
[25 



[26 

[27; 

[28 
[29 

[30 
[31 

[32 
[33 



M. Luke and M. J. Savage, Phys. Lett. B321 88 (1994). 

A. F. Falk, M. Luke and M. J. Savage, Phys. Rev. D53 6316 (1996). 

A. Kapustin and Z. Ligeti, Phys. Lett. B355 318 (1995). 

Z. Ligeti and Y. Nir, Phys. Rev. D49 4331 (1994). 

M. Gremm, A. Kapustin, Z. Ligeti and M. B. Wise, Phys. Rev. Lett. 77 20 (1996). 
V. Chernyak, Nucl. Phys. B457 96 (1995); preprint Budker INP-96-23, Eepj 



ph/9604376j 



M. Neubert, Phys. Rev. D46 1076 (1992). 

V. Eletsky and E. Shuryak, Phys. Lett. B276 191 (1992). 

P. Ball and V. M. Braun, Phys. Rev. D49 2472 (1994). 

I. Bigi, M. A. Shifman, N. Uraltsev and A. Vainshtein, Int. J. Mod. Phys. A9 2467 
(1994). 

M. Voloshin, preprint TPI-MINN-94/18-T. 

I. Bigi, A. G. Grozin, M. A. Shifman, N. Uraltsev and A. Vainshtein, Phys. Lett. 
B339 160 (1994). 

B. Barish et al, CLEO Collab., Phys. Rev. D51 1014 (1995). 



V. Gimenez, G. Martinelli and C. T. Sachrajda, preprint CERN-TH/96-175, |hcp- 
lat/9607055j . 



[34] E. Bagan, P. Ball, V. M. Braun e P. Gosdzinsky, Phys. Lett. B342 362 (1995) [E: 
ibid. B374 363 (1996)]; E. Bagan, P. Ball, B. Fiol e P. Gosdzinsky, Phys. Lett. B351 
546 (1995). 

[35] I. I. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. Lett. 71 
496 (1993). 

[36] M. Jezabek and J. H. Kuhn, Nucl. Phys. B320 20 (1989); A. Czarnecki and M. 
Jezabek, Nucl. Phys. B427 3 (1994). 

10 







m 







0] HI 




This paper 




0.1 


0.19 ±0.1 0.135 


0.44 ±0.25 


0.23 ±0.11 0.60 ±0.10 


-0.09 ±0.14 


0.66 



Table I. 

Results for the meson matrix element of the b quark kinetic energy operator. 
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Figure Captions 



Figure 1 

The 6-quark average square momentum K as a function of The straight line results 
from a two parameters fit. 

Figure 2 

The lepton spectrum in the semileptonic decay 5 — > X c £z/. 
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E lep t (GeV) 



Fig. 2 



